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We study the scaling behavior of the entanglement entropy of two dimensional conformal quantum 
critical systems, i.e. systems with scale invariant wave functions. They include two-dimensional 
generalized quantum dimer models on bipartite lattices and quantum loop models, as well as the 
quantum Lifshitz model and related gauge theories. We show that, under quite general conditions, 
the entanglement entropy of a large and simply connected sub-system of an infinite system with 
a smooth boundary has a universal finite contribution, as well as scale-invariant terms for special 
geometries. The universal finite contribution to the entanglement entropy is computable in terms 
of the properties of the conformal structure of the wave function of these quantum critical systems. 
The calculation of the universal term reduces to a problem in boundary conformal field theory. 
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I. INTRODUCTION 

The non-local correlations of a quantum mechanical 
system are encoded in the behavior of the entanglement 
properties of its wave functions. A pure quantum state 
of a bipartite system AU B defines a mixed state in the 
observed region A obtained from tracing out the degrees 
of freedom in the unobserved region B. The non-local 
correlations connecting regions A and B are encoded in 
the behavior of the von Neumann entanglement entropy, 
S = —TrpA In PA, where pA is the reduced density matrix 
of region A. The entanglement entropy of a local quan- 
tum field theory, relativistic or not is known to exhibit an 
"area law" scaling of the form S ^ p£^~^ in spatial di- 
mensions D > 1 where /x is a non-universal coefficient^'^. 
There has been growing interest in the scaling behavior 
of the entanglement entropy at quantum critical points 
and in topological phases. The entanglement entropy of 
quantum critical systems in _D > 1 should contain uni- 
versal subleading terms, whose structure for a general 
quantum critical system is not yet known. 

The scaling behavior of the entanglement entropy has 
only been studied in detail in quantum critical systems 
in _D = 1 space dimension. Such systems are described 
by a (1 -I- l)-dimensional conformal field theory (CFT). 
In a 1 -|- 1-dimensional CFT, the entanglement entropy 
of a subsystem A of linear size £ of an otherwise infinite 
system (i.e. of linear size L oo) obeys a logarithmic 
scaling law,— i^i^ii ~ f ^^ia)~^- ■ -j where c is the central 
charge of the CFT, and a is the short distance cutoff. 
There has been a number of studies on topics related to 
this ID logarithmic scaling form. For instance, a possible 
connection between this result and gravitational physics 
was suggested^. A similar logarithmic scaling behavior 
was found at infinite disorder fixed points of ID random 
spin chains^ii^. The quantum entanglement of quantum 
impurity systems has also been studied . -'^^i^^'^'^i-'^'^'^^ 

In this paper, we consider the universal scaling form 
of the entanglement entropy at 2D conformal quantum 



critical points (QCP) - two-dimensional quantum criti- 
cal systems with scale-invariant many body wave func- 
tions. At a 2D conformal QCP, equal-time correlators 
of local operators coincide with the correlation func- 
tions of an appropriate 2D classical system at criticality 
(which is described by an Euclidean 2D CFT)!^. The 
entanglement entropy of 2D conformal QCPs was first 
considered in Ref.[17[, where a scaling form was found: 
S — p,l — |(Ax)ln(£/a) -I- . . ., where c is the central 
charge of the 2D Euclidean CFT associated with the 
norm squared of the wave function and Ax is the change 
of the Euler characteristic x. Ax — Xaub~Xa^Xb- No- 
tice for a region A <Z B with a smooth boundary. Ax — 
and hence the logarithmic term vanishes. Hence, if region 
A has a smooth boundary, there is no universal logarith- 
mic term. In this case, we will show that instead there is 
a finite, 0(1), universal term jqcp in the entanglement 
entropy at these quantum critical points, i.e. 



QCP 



IJ-i + IQCP + 



(1.1) 



Through explicit calculations and using general argu- 
ments based on CFT, we will show that Jqcp has a 
topological meaning in the sense that it is determined 
by the contributions of the winding modes of the under- 
lying CFT. 

In a topological phase in 2D, the entanglement entropy 
scales aai^ii^ 



topo 



'ytopo 



0{£-'), 



(1.2) 



where q; is cL non-universs-l coefficient cind ^topo 5 

the 

topological entanglement entropy, is a topological invari- 
ant, the logarithm of the so-called total quantum dimen- 
sion V of the underlying topological field theory describ- 
ing the topological phaseJ^ii^ Topological phases have 
non-trivial ground state degeneracies on surfaces of non- 
trivial topology. The topological entanglement entropy 
7topo also depends on the global topology of the man- 
ifold, and on surfaces with non-trivial topology, on the 
degenerate ground state on that surfacci^ 
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Although superficially similar, the finite universal con- 
tributions to the entanglement entropy in topological 
phases and conformal quantum critical points, 7topo and 
IQCP, have a different origin and structure. In the case of 
a topological phase, 7topo is in general determined by the 
modular 5-matrix of the topological field theory of the 
topological phase fiSiii^i^ This modular 5'-matrix governs 
the transformation properties of the (degenerate) ground 
states of the topological phase on a torus under modular 
transformations, r — 1/t, where r is the modular pa- 
rameter of the torus.— However, we show below that for a 
general conformal quantum critical point, whose ground 
state wave function is given by the Gibbs weights of a Eu- 
clidean rational unitary CFT, the universal term ^qcp 
is determined by the modular S'-matrix associated with 
the norm squared of the wave function. Thus, the mod- 
ular S'-matrix of the topological phase and that of the 
wave functions of 2D conformal quantum critical points 
have a conceptually different origin. In particular, in all 
the cases we checked here, Jqcp and 7topo contribute 
with opposite signs to their respective entanglement en- 
tropies, as implied by the conventions we used in Eq. lfTT]) 
and Eq. fL^ . 

We will show that, when the logarithmic terms in the 
entanglement entropy cancel, the finite terms ^qcp are 
universal and are determined not only by the central 
charge but also by the restrictions on the states im- 
posed by the compactification conditions. Furthermore, 
the form of the result for the entanglement entropy of 
Ea. (|2.6[) implies a connection with boundary CFT, as 
developed by Cardy .^^'^'^ Thus, in addition of it being 
determined by the central charge c, it must also depend 
on the operator content of the CFT. For the same reason, 
the structure of Ea. l|2.6p also suggests a direct connection 
between this problem and the Affleck-Ludwig boundary 
entropy of ID quantum CFTs.-^ 

The paper is organized as follows. In Section [ill we ap- 
ply this approach first to the simpler case of the quantum 
Lifshitz model (and the related quantum dimer models, 
QDMs) on planar, cylindrical and toroidal geometries. 
These results apply to the QCPs of (generalized) quan- 
tum dimer model on bipartite lattice o^^'^^i^'^'^^i^^'^" and 
in quantum eight-vertex models^. Through explicit cal- 
culations for various geometries, we show that that, when 
the logarithmic terms in the entanglement entropy can- 
cel, and that the subleading finite terms jqcp are uni- 
versal, determined not only by the central charge but 
also by the restrictions imposed by the compactification 
conditions. In Section lllll we generalize this result to 
all 2D conformal QCPs whose scale-invariant wave func- 
tions have norms that are the partition functions of 2D 
Euclidean Rational CFTs (RCFT), CFTs with a finite 
number of primary fields^"^. More specifically, we show 
that the finite term in the entanglement entropy of the 2D 
wave function is determined by the change of the Affleck- 
Ludwig boundary entropy of the ID CFT - a quantity 
determined by the modular S'-matrix of the associated 
CFT and by the coefficients in the fusion rules. We also 



discuss specific examples of this class including 2D quan- 
tum loop models^ which, with the naive inner product, 
are known to be quantum critical . ^"^'^^ We also briefly 
discess the quantum net modelsj^i^'^i^^ In Section IIVI 
we conclude with a summary and a discussion on open 
questions. In particular, we comment on the implications 
of our results to the nature of related topological phases. 



II. QUANTUM LIFSHITZ MODEL 
UNIVERSALITY CLASS 

The quantum Lifshitz model^S^ (QLM) in two space 
dimensions is defined by the following Hamiltonian with 
an arbitrary paramter k: 



H = d^x 



n' 



1 



(2.1) 



where (pis a, scalar field T\ ~ (pis its canonical momentum 
conjugate to 0. The QLM Hamiltonian Eq. (|2.ip defines 
a class of QCP's with dynamic critical exponent z — 2, 
and a continuous parameter k. 

This remarkable property of the model is evident in the 
exactly known wave function for the ground state ICS') 
which is a superposition of all field configurations (/)(x, y) 
with the configuration dependent weig 



■^Gs[4>] = m\GS) 



with 



S[0] 



A (v0(.) 



and the norm squared of the state 



Z=\\^ 



GS\ 



= D(j) e 



(2.2) 



(2.3) 



(2.4) 



Notice Z is identical to the partition function for the 
Gaussian model, which defines free boson Euclidean 
CFT^^, albeit with the "stiffness" k. Hence Eq. (|2l|) 
defines an infinite class of 2D conformal QCP's all as- 
sociated with free boson CFTs. 

The QLM can be viewed low energy effective field the- 
ory capturing universal aspects of various microscopic 
lattice models with (p playing the role of coarse grained 
height fieldi^i^i^S with the "stiffness" k determined by 
the appropriate "microscopic" coupling constant s^^'^^ . 
For such a mapping to work, the constraints of the lat- 
tice models should be build in through compactification 
of the boson field by demanding all physical operators 
to be invariant under the shift oi cj) ^ (p ~\- 27rr or equiv- 
alently all physical operators to take the form of vertex 
operators e™'^/'' for integer n. In subsection III CI we will 
discuss specific examples of this mapping corresponding 
to particular values of k using the convention of fixing 
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r = 1. The examples will include so-called Rokhsar- 
Kivelson point (RK) of the quantum dimer modeP^ and 
its generalizations^fliiSiiil and the quantum eight-vertex 
model^ special choices of the Baxter weight*^. Since 
k can be varied in the QLM, this theory has an exactly 
marginal operator, resulting in continuously varying criti- 
cal exponents (scaling dimensions) of the allowed (vertex) 
operators 4i 



A. Entanglement entropy and partition functions 
for 2D conformal QCPs 

To investigate the universal finite terms in the entan- 
glement entropy at 2D conformal QCPs, we will rely 
on the approach described in the work of Fradkin and 
Moore.— They showed that trp^, where pA is the (nor- 
malized) reduced density matrix of a region A, with 
A C B separated by the boundary T, for the ground 
state 4'o on A(J B, is given by 



take the form 



= — 



/ rr rj \ n—1 



(2.5) 



Here Zn is the partition function of n copies of the equiva- 
lent 2D classical statistical mechanical system satisfying 
the constraint that their degrees of freedom are identi- 
fied on the boundary F, and is the partition function 
for n decoupled systems. The partition functions on the 
r.h.s of Eq. (|2.5p are Za — 11^*0 IP ^i^h support on region 
A and with support in region B, both satisfying 

generalized Dirichlct {i.e. fixed) boundary conditions on 
F of A and _B, and Za\jb = H^olP is the norm squared 
for the full system. The entanglement entropy S is then 
obtained by an analytic continuation in n, 

S = -tr (pa In pa) 
= - lim ^trp4 



= -log 



ZaZp 



^AUB 



(2.6) 



Hence, the computation of the entanglement entropy is 
reduced to the computation of a ratio of partition func- 
tions in a 2D classical statistical mechanical problem, an 
Euclidean CFT in the case of a critical wave function, 
each satisfying specific boundary conditions. 

In order to construct trp^, we need an expression 
for the matrix elements of the reduced density ma- 
trix {(t)'^\pA\(t)'^) ■ Since the ground state wave function 



Eqs. (|2.2p and (|2.3p is a local function of the field (/)(a:), 
a general matrix element of the reduced density matrix 
is a trace of the density matrix of the pure state ^'gs'[0] 
over the degrees of freedom of the "unobserved" region 
B, denoted by (j)^{x). Hence the matrix elements of pA 



(2.7) 



where the degrees of freedom satisfy the boundary condi- 
tion at the common boundary F: 



BCr 



ir = 



ir = 



(2.8) 



Proceeding with the computation of trp^ , it is immediate 
to see that the matrix product requires the condition 



for i 



1, 



and 



trace condition. Hence, trp^ takes the form 



from the 



^Wa = ^ 



1 



-ill XlD^tD<l>l 



e~S:?=i(s(0f)+s(0f)) 
(2.9) 



subject to the boundary condition BCr of Eq. (|2.8p . No- 
tice that the numerator, Z„ is the partition function on 
n systems whose degrees of freedom are identified in F 
but are otherwise independent. Also notice the absence 
of the factors of 1/2 in the exponentials of Eq. (|2.9p . 

The other important consideration is that the com- 
pactification condition requires that two fields that differ 
by 2-Kr be equivalent. Hence, the boundary condition of 
Ea. (j2.8p is defined modulo 2i:r. (Equivalently, the proper 
form of the degrees of freedom is e''^.) This means that 
one can alternatively define Z„ as a partition function 
for n systems which are decoupled in the hulk but have 
a boundary coupling of the form (in the limit Ar — > oo, 
which enforces the boundary condition) 



Ar ^ cos( 

i=l 



(2.10) 



Here the fields (f)i extend over the entire region AU B. 
Thus, this problem maps onto a boundary CFT for a sys- 
tem with n "replicas" coupled only through the boundary 
condition on the closed contour F, the boundary between 
the A and B regions. 

For the special case of the free scalar field, one can 
simplify this further by taking linear combinations of the 
replica fields. Then the condition that the scalar fields 

agree with each other on F can be satisfied by forming 
n—1 relative coordinates tpi = (f>i — (i — 1, . . . ,n—l) 
that vanish [mod 2'Kr) on F, and one "center of mass 
coordinate" field 4> = X]"=i 4'i ^^sX is unaffected by 
the boundary F (reflecting the fact that nothing physi- 
cal takes place at F). Hence, the computation of trp^ 
reduces to the product of two partition functions: 
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1. The partition function for the "center of mass" field 
4>] since 4> does not see the boundary F, this is just 
the partition function Zaub for a single field in the 
entire system. 

2. The partition function for the n — 1 fields ipi which 
are independent from each other and vanish {mod 

2'iTr on r. We denote this by (^r*)" • However, 
the fields (pi on the A and B regions are effec- 
tively decoupled from each other. Hence, this par- 



tition function further factorizes to 



rjD yD 



where Z^ and Z^ are the partition functions for 
a single field on A and B respectively, satisfying 
in each case Dirichlet (fixed) boundary conditions 
{mod 2'Kr) at their common boundary T . 

Thus, we can write the trace trp^ as 



Za 



UB 



Z 



AUB 



zl 

Zf 



yD yD 
^A ^B 



Zaub 

(2.11) 

Here the denominator factor, Z''^^^^ comes from the nor- 
malization factors, and represents the partition function 
over the entire system. The entanglement entropy is 
ther 



,17 



B. The Entanglement Entropy of the Quantum 
Lifshitz Universahty Class 

Here we calculate jqcp at QCPs of the QLM univer- 
sality class defined by Ea. l|2.ip for three different geome- 
tries: (i) a cylindrical geometry, (ii) a toroidal geometry, 
and (iii) a disk geometry. For the cylinder and disk we 
assume the Dirichlet boundary conditions at the open 
ends. We use the known results on the free boson parti- 
tion function (j2.4p for different topologies and boundary 
conditions'^ii^^i^i^i^i^, which are necessary for the cal- 
culation of entanglement entropy. It is useful to note that 
the action Eq. (|2.3p for general value of the "stiffness" k 
turns into the standard form: 



1 

8^ 



d^x {dfj^ifiY 



(2.15) 



upon a rescaling of the field w2k(f) — '-p. If is com- 
pactified with radius r = 1, the rescaled field iy9 has an 
effective compactification radius R — v2kr^. We find 
iQCP to depend linearly on In R in all cases we consider. 



1. The Cylinder 



log Z^ -log Z^ + log Zaub 



Faub, 
(2.12) 

which, as indicated in the r.h.s of Eq. (|2.12p reduces to 
the computation of the free energies , and Favjb, 
for the equivalent 2D Euclidean CFT on regions A and B, 
each satisfying Dirichlet (fixed) boundary conditions on 
the common boundary F, and on the full system, A U B, 
respectively. 

The behavior of the free energy of a CFT as a func- 
tion of the system size I has been studied in detail. The 
divergent terms, as i ^ co, have the form^i^i^ 



F{1)^ M^+at-lx^n{- 
6 V a 



0{l) 



(2.13) 



provided the boundary F is smooth (and differentiable). 
Here, /o and a are two non-universal quantities, and a is 
the short-distance cutoff; c and x are, respectively, the 
central charge of the CFT and the Euler characteristic of 
the manifold. It follows from this result that the entan- 
glement entropy for region A takes the fornJ^ 



6 V a 



0{l). 



(2.14) 



provided the boundary F is smooth. In all the geometries 
we discuss, the change in the Euler characteristic van- 
ishes. Ax = 0, and there is no logarithmic term. How- 
ever we will show below that, if the logarithmic terms 
cancel, there exist a universal finite 0{\) term, as well 
as other universal dependences on the geometry (such as 
aspect ratios) . We will now extract these universal finite 
terms. 



Dirichlet 

V 



Dirichlet 

r 




->•-<- 



L, 



FIG. 1: Cylinder 

Let us begin by considering first a system on a long 
cylinder of linear size L and circumference ^ with L 3> 
i. Region A to be observed, is a cylinder of length La 
and circumference The complement region, B, is a 
cylinder of length Lq (see FigH]), also with circumference 
t. We assume that the QLM wave function Ea. (|2.2p and 
hence the associated 2D partition function Eq. (|2.4p obey 
the Dirichlet boundary conditions at both ends of the 
cylinder, A U B. 

From Eq. (|2.12p . the entanglement entropy Sa — Sb ^ 
S is given by 

s = -inz^ij{LA,e)-iiizEDiLB,e)+iiiz^'if{LA+LB,e) 

(2.16) 

Here Zdd{L,£) is the partition function of Eq. ()2.4p for 
a boson with compactification radius R on cylinder of 
length L and circumference £ with Dirichlet boundary 
conditions on both ends, which is well known 4^ 



ZDD{L,e)^u 



R rj{q^) 



(2.17) 
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where R = V2r^k is the effective compactification ra- 
dius (as before), and A/" is a non-universal regularization- 
dependent prefactor, responsible for the area and perime- 
ter dependent terms in the free energy shown in 
Ea. (|2.13p . (There are no logarithmic terms for a cylin- 
der or a torus as their Euler characteristic x vanishes.) 
In Ea. (j2.17p r = ij is the modular parameter, encoding 
the geometry of the cylinder, and q — e^'^*'^. The elliptic 
theta-function ^3{t) and the Dedekind cta-function ri{q) 
are given by 

oo oo 
n— — 00 n— 1 

The important feature of Eq. (j2.17|) is the factor 1/R, 
the contribution of the winding modes of the compacti- 
fied boson on the cylinder with Dirichlet boundary con- 
ditions. 

Putting it all together, it is straightforward to find an 
expression for the entanglement entropy using Eq. (j2.6p . 
In general, the entanglement entropy depends on the ge- 
ometry {e.g. the aspect ratios L/£) of the cylinders, en- 
coded in ratios of theta and eta functions. However, in 
the limit La ^ (, in which the length of the cylinders are 
long compared to their circumference, the entanglement 
entropy given by Eq. (|2.16p and Ea. (|2.17p takes a simple 
form 

S = ^ie + \nR, (2.19) 

where /i is a non-universal constant that depending on 
the regularization-dcpcndcnt pre-factor J\f of Eq. (|2.17p . 
Hence, there is a 0(1) universal contribution to the en- 
tanglement entropy jqc p = In i? for the cylinderical ge- 
ometry. The explicit dependence of Jqcp on the effective 
effective compactification radius R = V2fcr^ shows that 
it is determined by the winding modes of the compacti- 
fied boson and thus it is a universal quantity determined 
by the topology of the surface. In particular we find that 
the universal piece of the entanglement entropy, ^qcp, 
for a compactified boson is a continuous function of the 
radius R, a consequence of the existence of an exactly 
marginal operator at this QCP. We find the similar rela- 
tions for all topologies we considered. We will come back 
to this point in section III C[ in the context of several 
microscopic models of interest. 

2. The Torus 

We now consider the case in which the full system 
AU B is a torus for which the real part of the modu- 
lus L/£ ^ 1, as shown in FigO The two subsystems, 
A and B are now two cylinders, of length La and Lb 
respectively {L — La + Lb), both with the same circum- 
ference i. We will thus need the partition function on a 
torus and on two cylinders (with both ends of the cylin- 
ders obeying Dirichlet boundary conditions.) The trace 



trp^ now becomes 

The partition functions for the two cylinders, A and B 
has the form of Eq. (I2.17p . The partition function for 
the torus i s-"'"" 

^toius(i,^) = ^^j^lhidor ^"^7^^^ > (2-21) 

where Z^^^^^^{^,£) is the partition function on a cylin- 
der of length -J and circumference i, with Neumann 
boundary conditions at both ends: 

«.,.(f,^)^A^;?^, (-^) 

where r = i j and q = exp(27rir). 

In the limit La ^ ( ^ a and Lb ^ i ^ a, the 
entanglement entropy for the toroidal geometry is 

S = ^ii + 2ln (^^^ . (2.23) 

Hence, for the toroidal geometry, the universal term is 
7QCP = 21n(fcr2) = 21n(i?V2). In Eq.jMS]) is, once 
again, a non-universal factor which depends on both the 
short distance regularization and boundary conditions (in 
fact, it is not equal to the constant we also called in 
the entanglement entropy for the case of the cylinder, 
Ea. (|2.19p .) As was the case for the cylindrical geometry, 
in the case of the torus "fqcp is also determined by the 
contribution of the zero modes of the compactified boson 
to the partition fimctions. Thus, here too, jqcp depends 
on the effective boson radius R = v2fcr^. However, the 
different values of jqc p in Eq. (|2.23p and Eq. (|2.19p is due 
to the fact that on the torus all three partition functions 
have contributions from the zero modes. 




FIG. 2: Torus 
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The Disk 



tion functions are 




FIG. 3: Disk 

Finally, we compute the entanglement entropy for the 
disk geometry, shown in FiglS] The line of argument used 
above applies here as well. This is the case discussed in 
Ref. [l3| , where it was found that the logarithmic term in 
the entanglement entropy cancels exactly if the boundary 
r is smooth. Here we compute the (subleading) finite 
universal piece. 

To compute the entanglement entropy we need to com- 
pute three partition functions, on the two disks A and 
AUB, and on the annulus B, all with Dirichlet boundary 
conditions. These partition functions were computed in 
the literature long ago for an uncompactified boson4^i^ 
They can be obtained from the partition functions on 
cylinders, with Dirichlet-Dirichlet (for the annulus) and 
Dirichlet-Neumann (for the disks) boundary conditions 
by a conformal mapping w = -^Inz, from the z complex 
plane to the cylinder (labeled by w). The partition func- 
tion for the annulus (region B) of inner circumference 
£ and outer circumference L (with Dirichlet boundary 
conditions) is 



1 ^3(^) 



(2.24) 



Except for the factor of l/v2fcr^, which is due to 
the zero modes of the compactified boson, this result 
agrees with those of Ref. p7| . In Ea. l|2.24p we have 
used qs = e^'^"'^ = j- (with the modular parameter 
ln(f)). 



2tt 



Similarly, the partition functions on the two disks, re- 
gions A and A U B, are conformally mapped to two in- 
finitely long cylinders (as the UV cutoff a — > 0) with 
Neumann-Dirichlet boundary conditions. These parti- 



idisk 



5/12„l/4 ^4 (t) 
TT 



(2.25) 



where q = (f)^:(f:)^ for regions A and A U B, respec- 
tively, and T is their corresponding modular parameter; 
•diij) is the elliptic theta- function 



(2.26) 



The resulting entanglement entropy for the planar 
(disk) geometry is found to be 



S= - \n 
2 



Hi 



Ini?. 



(2.27) 



Hence, for the case of the disk there is also a univer- 
sal finite piece in the entanglement entropy, jqcp = 
In y/2kr^ = InR. As in the cases discussed above (the 
cylinder and the torus), here too ^qcp has a topologi- 
cal origin as it is due to the winding modes of the com- 
pactified boson. However, unlike the case of the of the 
cylinder and toroidal geometries, in the case of the disk 
there is also a dependence on the aspect ratio L/i (th e 
double logarithmic term), as already noted in Ref.[l7|. 
(Note that we included the factor of I/tt in the double 
logarithm since it arises from the conformal mapping.) 



C. Entanglement Entropy of Quantum Dimer 
Models and Related Systems 

The results on the entanglement entropy of the pre- 
ceding subsections apply to several "microscopic" sys- 
tems of interest. The simplest of them is the quantum 
dimer model on bipartite lattices at the RK point (as- 
sociated with the RK wave function of the QDM). As 
noted in Ref.[l6|, the RK point of the QDM maps onto 
the quantum Lifshitz model for a particular value of the 
radius r = 1 and stiffness fc = 2 (in the notation used 
here.) This corresponds to a 2D Euclidean boson CFT 
at the free fermion radius. Of course, this is not an acci- 
dent, since in this case the lattice partition functions can 
also be computed exactly by pfaffian methods,— i^i^ and 
hence it is a free Dirac fermion system. 

Generalized quantum dimer models have been dis- 
cussed recently .^^'^°i"^° In these models the wave func- 
tions correspond to dimer models with weights that de- 
pend on the number of dimer pairs on the plaquettes. 
For a considerable range of values of these weights the 
system remains critical and can also be mapped onto a 
quantum Lifshitz model, albeit with a different stiffness 
connected with the presence of an exactly marginal op- 
erator. Thus, in these models the stiffness varies con- 
tinuously as a function of the microscopic weights. This 
dependence, discussed in detail in Ref.jl^, is of course 



7 



non-universal, as it depends on the microscopic structure 
of the system. Nevertheless, the critical exponents have a 
universal dependence on the stiffness. The same applies 
to the universal piece of the entanglement entropy jqcp, 
which can be read-off from the results presented in this 
section. 

Similarly, the quantum eight-vertex model wave 
function^^ also maps onto a free fermion problem for a 
special choice of weights4^ For general values of k the 
fermions are interacting (see the discussion below) but 
the effects only enter through an exactly marginal opera- 
tor. The mapping of the quantum 2D eight-vertex model 
to the quantum Lifshitz model was shown in detail in 
Ref. [iBl where the relation between the stiffness k of the 
compactified boson and the Baxter weights is given ex- 
plicitly, k and the weight c in the Baxter wave function 
(along the six vertex line) are related by 



TT 

2k 



cot 



1 



(2.28) 



for a boson with compactification radius r ~ 1 or, equiv- 
alent, an effective radius R = v2kr^. 

The results of the preceding subsections on the entan- 
glement entropy for the quantum Lifshitz model apply 
to the lattice models almost without change. Once the 
mapping of the stiffness to the microscopic parameters 
(as in the case of the quantum eight vertex model) is 
known, the universal piece, jqcp, can be read-off imme- 
diately. The only caveat here is that in lattice models it is 
impossible to have closed simply connected regions with 
smooth boundaries. The resulting paths of the effective 
coarse grained quantum Lifshitz model will always have 
singularities, such as corners, which contribute with a 
logarithmic dependence to the entanglement entropy (as 
discussed in Ref. [l7| ) rendering the finite terms gener- 
ally non-universal. The cylinder and torus geometries 
are exceptional in this sense, and allow for a direct check 
of these ideas in microscopic models, either through an 
exact solution or by means of numerical computations. 

We end this discussion by giving the results for the 
universal entanglement entropies ^qcp for the Lifshitz 
universality class at the free fermion (or dimer) and 
Kosterlitz-Thouless transition of the dimer and Baxter 
(six vertex) wave functions for all three geometries. (See 
the summary of Table [H) At the "free dimer" point 
(the free fermion point of the dimer models) the stiffness 
k — 2 (corresponding to = \/2 in the Baxter wave 
function), and the universal term of the entanglement 
entropy for a disk geometry is 7qcp ~ lnv2fcr^ = ln2. 
For the cylinder, also at the free dimer point, we also 
found 7qcp^°'^ = In 2, while for the torus we obtained 
7qpp = 2 In 2. (Below we will discuss the relation of 
these results with the topological entanglement entropy 
of the nearby Z2 topological phase.) Away from the free 
dimer (or fermion) points, the stiffness k changes and so 
does the entanglement entropy. Thus, at the Kosterlitz- 
Thouless transition point of both the dimer and six ver- 
tex wave functions (where the Baxter weig ht is c = \/2), 



the stiffness is fc = 1. (At this point the associated c = 1 
CFT has an SU (2)i Kac-Moody current algebra, and the 
effective compactification radius here is i? = \/2-) The 
(finite) entanglement entropies now are 7qcp = 2 In \/2, 
l';^^f''' = 0, and7^^'V-lnN/2. 



R 


cylinder 


torus 


disk 


2 (RKpoint) 


In 2 


2 In 2 


In 2 


^2 (KTpoint) 


Inv^ 





In V2 



TABLE I: Universal entanglement entropies 'yqcp of the lat- 
tice models in QLM universality class in the cylinder, torus, 
and disk geometries, yqcp based on calculations from QLM 
is quoted at the free fermion point (or RK point) R = 2, and 
at the Kosterlitz-Thouless {SU{2)i) point, R ^ V2. 

The only caveat in applying the calculation of 'Jqcp 
in the QLM to microscopic models is that is impossi- 
ble to have closed simply connected regions with smooth 
boundaries on a lattice. Hence the resulting paths of 
the effective coarse grained QLM will always have singu- 
larities (such as corners) which contribute a finite loga- 
rithmic dependence to the entanglement entropy.— The 
cylinder and torus geometries are exceptional in this 
sense, and allow for a direct check of these ideas in mi- 
croscopic models, either through an exact solution or by 
means of numerical computations. 



GENERALIZED CONFORMAL QCPS 
ASSOCIATED WITH RCFT 



III. 



We now generalize the application of Eq. (|2.6p to the 
computation of the entanglement entropy to more gen- 
eral case of conformal QCPs, specifically those associated 
whose wave functions have an associated 2D Euclidean 
RCFT (a CFT with a finite number of primary fields.) 



A. Entanglement entropy and Boundary 
Conformal Field theory 

The ground state wave function for a conformal quan- 
tum critical point can be expressed as Gibbs weight as- 
sociated with a 2D Euclidean CFT: 



1 



-5[0]/2 



(3.1) 



as in the case of the QLM discussed in the previous sec- 
tion. Hence there is a one-to-one mapping between the 
norm square of the wave function and the partition func- 
tion of a local 2D Euclidean CFT, and also between the 
equal-time correlators of the operators of the 2D confor- 
mal QCP map onto and the correlators of primary fields 
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of the 2D Euclidean CFT. Furthermore, we will also as- 
sume that the associated Euclidean CFT is unitary (the 
^-matrix to be defined below is unitary) and that it is 
a RCFT. The restriction to unitary RCFT allows us to 
exploit well developed technology for this large class of 
CFTs^i^^, especially that of operator product expansion 
(OPE) and of modular S-martirx, in calculation of jqcp- 
The behavior of RCFTs with specified boundary condi- 
tions (especially their partition functions), is the subject 
of boundary conformal field theory, and was discussed 
extensively by Card y^^i^^ . We will follow the approach 
and results of Cardy in this section. We also need to 
specify the boundary conditions at the ends of the cylin- 
der, i.e. the boundary states of the boundary CFTi22 
Let us denote these conformal boundary conditions by 
(a,/3). The associated (conformally invariant) bound- 
ary states (a| and \b) can be constructed for each CFT. 
On the other hand, at the common boundary F between 
the regions A and B, all ri — 1 fields must obey fixed 
('Dirichlet') boundary conditions. As shown by Cardy^^^ 
this boundary condition is quite generally given by the 
boundary state |0) in the conformal block of the identity 
1. 

For simplicity, we will consider here only the geome- 
tries of a cylinder (with specific boundary conditions at 
each end) and a torus. As in Ea. (|2.6p we will need to 
compute the free energies of region A, B and AyjB with 
fixed boundary conditions. 

The partition function for a RCFT on a cylinder of 
length L and circumference i?, with boundary conditions 
a and b on the left and right ends respectively, Z^/b, can 
be expressed in terms of the characters Xi of the RCFT: 



Za/b — 



Ni 



abXj 



-Tr£/L^ 



(3.2) 



where the integers iV;^^ are the fusion constants, the co- 
efhcients in the OPE of the RCFT, 



$„ X $b = > ' NL^ 



3 



(3.3) 



The Virasoro characters Xj s-re given by the trace over 
the descendants 1$^) of the highest weight state, which 
are obtained by acting on it with the Virasoro generators 
L_„ {n > 0): 



7r£/L\ _ iTlc/2'iL 



X,(e— /-) = e 



tra (e^^^°) : (3-4) 



where c is the central charge of the CFT, Lq is the 
n = Virasoro generator. Here the modular parameter is 
T = iijlh. Under a modular transformation r —XjT^ 
which exchanges the Euclidean "space" and "time" di- 
mensions of the cylinder (i.e. it flips the cylinder from 
the "horizontal" to the "vertical" position), the charac- 
ters transform as 



.11 L\ ^ S3 



(3.5) 



where Sf is the modular S-matrix of the RCFT. The 
modular 5'-matrix and the fusion coefficients are related 
by the Verlinde formula— 



ci ci cb 
ab / J ci 



(3.6) 



The limit of interest here is, once again, L ^ £. Un- 
der a modular transformation, the partition function of 
Eq. p.2p becomes 



Za/b = Y.^abSiXj{e 



-47ri/A 



(3.7) 



In the limit j- ^ 0, Za/b is dominated by the the de- 
scendants of the identity 1 (up to exponentially small 
corrections). Hence, in this limit. 



Za/,^Y.NabS?Xo {e-^'^'") 



£67 



and In Za/b becomes 



\nZ, 



a/b 



ttLc 



(3.8) 
(3.9) 



dropping UV singular (non-universal) terms. The quan- 
tity In gab in Eq. (|3.9p is the boundary entropy of a bound- 
ary RCFT introduced by Affleck and Ludwig— , where 
the "ground state degeneracy" gab is given by 



(3.10) 



Using Eq. (|2.6p . these standard results imply that the 
entanglement entropy of the 2D rational conformal QCP 
for a cylindrical geometry(see Fig[T]). For boundary con- 
ditions a and b at the two ends associated with regions 
A and B, the entanglement entropy is 



5* 



In 



yaO yOb 
yab 



Ai^-ln 



Y.iKbS'l 



V 9ab 



(3.11) 



where we explicitly used the fact that the state at the 
common boundary F should be fixed to be the fixed BC 
with boundary state |0). 

The result Eq. ()3.11|l provides an explicit way to com- 
pute ^Qcp for the entire class of many-body wave func- 
tions at QCPs associated with RCFT in terms of the data 
of the RCFT: 



IQCP 



In 



El Kb s". 



(3.12) 
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This is the main result of this section. It shows that 
IQCP is in general determined by the OPE coefficients 
(which encode the boundary conditions on the parti- 
tion functions) and by the modular S'-matrix, Si, of the 
RCFT associated with the norm squared of the many- 
body wave function at the given QCP. 

It is important to note that it is also possible to de- 
fine a unitary S'-matrix that governs the transformation 
properties of the wave function itself under a modular 
transformation. This modular S'-matrix plays a central 
role in 2D topological phases and in topological field 
theories ji^i^ii^ However, only for topological theories 
these are two S-matrices are the same and in general 
they are different or not even defined at all. We will 
come back to this issue in the discussion section. 

A particularly simple result is obtained for the case of 
a cylinder with fixed boundary conditions on both ends. 
In this case, Za, Zb and Zaub are cylinders with fixed 
boundary conditions, and hence the boundary states for 
all three cases are in the conformal block of the identity 
1. Since in this case the only non- vanishing OPE coeffi- 
cient is A^QQ = 1, the universal term of the entanglement 
entropy, jqcp, depends only on the element Sq of the 
modular S-matrix of the RCFT: 



IQCP 



lnS°. 



(3.13) 



For the case in which the full regionylU-B is a torus, we 
can use an analogue of Eq. (|3.1ip by writing the partition 
function Zaub in the denominator of Eq. (j3.1ip as a mod- 
ular invariant. In the limit of interest L ^ the denom- 
inator gab of Eq. (|3.1ip is replaced by a sum of terms with 
similar structure corresponding to a sum over boundary 
conditions (and twists) needed to represent the torus (see, 
for instance, Ref.jl^). Similarly, Eq. (|3.11[) can also be 
applied to the disk geometry upon a conformal mapping 
as it was done for the case of the compactified boson in 
section III B 31 



B. Applications 

We will now discuss some examples of interest. In 
applying the results Ea. (|3.12p to specific systems, one 
should keep in mind that that choice of the inner prod- 
uct of the 2D quantum theory can play a subtle role. As 
it was pointed out recently by Fendley— , a scale invariant 
wave function does not necessarily imply scale invariance 
of the correlators. Their actual behavior depends also 
on the choice of inner product. Here we have assumed 
that the states labeled by the set of field configurations 
4>{x, y) form an orthogonal basis. Hence, the norm of the 
wave function is a sum over states with the local weights 
squared. However what matters is that the matrix ele- 
ments (and in particular the norm of the states) be scale- 
invariant. A number of interesting counterexamples are 
knowni^ The QLM is a special case where such "naive" 
inner product maintains scale invariance. This is due to 
the existence of exactly marginal operators in the QLM. 



Below we discuss four cases where the ground state 
wave function with the "naive" inner product describes 
QCPs: (i) a QCP associated with the 2D Ising CFT, (n) 
the QCPs associated with compactified boson CFT, (iii) 
QCPs in quantum loop models^i^, and (iv) quantum 
net models^i26,56^_ (ggg footnote Ref.jsl].) 



1. The 2D Ising wave function 

As an example of a system described by an RCFT we 
consider a 2D quantum spin system whose ground state 
wave function has for amplitudes the Gibbs weights of 
the 2D classical Ising model. This system is quantum 
critical if the square of the weights (which also have the 
form of a Gibbs weight for the 3D Ising model) are at the 
critical point of the 2D Ising model, the Onsager value. 

The critical point of the 2D Ising model is the simplest 
RCFT. It has central charge c — 1/2, and three (bulk) 
primary fields: 1) the identity (1, with conformal weight 
h — Q), 2) the energy density (e, with conformal weight 
h = 1/2), and 3) the spin field (cr, with conformal weight 
1/16), which obey the operator algebra (OPE) 

e X e = 1 
e X cr = (7 

axa^l+£. (3.14) 

The critical Ising model has three possible boundary 
states:^'^ 1) the spin up state |-|-), 2) the spin down state 
|— ), and 3) the free state |/). (Either the up or the 
down state can be regarded as the fixed boundary state.) 
These three boundary states, |-|-), |— ), and |/) are in the 
conformal blocks of the identity 1 (denoted by |0)), the 
energy density e (denoted by and the spin field a 

(denoted by respectively. The boundary states are 

given hy^ 



I/) ^ 



The modular S-matrix is 

s=| I 

\ 72 



|0) = 










^> 












'l6 


= |o)- 


|£>- 





I 
1 
\_ 

2 

1 



(3.15) 



(3.16) 



where the columns are labeled by the highest weights 0, 
1/2, and 1/16, in that order. 

The entanglement entropy for this wave function can 
now be computed, using the result of Eq. (|3.1ip . We will 
take region A U S to be a long cylinder of length L and 
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circumference and regions A and B to be two cylin- 
ders of lengths and Lb respectively, with the same 
circumference £, and with L = La + Lb- 

Let us take the boundary conditions at both ends of 
A U B to be free. By a conformal mapping, this maps 
onto the disk. Back on the cylinder, the free boundary 
condition is described by the boundary state |/), which 
is in the conformal block of the primary field a. On the 
other hand, at the boundary F between regions A and B, 
we have the fixed boundary condition, the up state 1+) 
We readily find 



5<T,o — N^^qS" 



1 



5.,.-<.5o" + iV^,^5^? = l. (3.17) 
The universal term of the entanglement entropy, ^qcp 



now IS 



IQCP 



= - In 



9a090b 
9ab 



In 



S2 



In 2. 



(3.18) 



On the other hand, we could consider instead the case 
of fixed boundary conditions at both ends of the cylinder 
AUB. This corresponds to the boundary state |0). Since 
the boundary condition on T is always fixed, "/qcp is now 



IQCP = - In 5*5 = In 2. 



(3.19) 



In the case where AUB is torus of large circumference L 
and small circumference £ (hence with modular parame- 
ter T = i£/L), the regions A and B are cylinders each of 
length La and L b and circumference with fixed bound- 
ary conditions at both ends. The partition function for 



the torus, Z^"^", 



^AUB 





Mr) 




Mr) 




Mr) 








+ 


V{r) 


+ 


tj{t) 


) 



(3.20) 



Using the modular invariance of Z on the torus (r — > 
-1/r), one finds that in the hmit L > ^, Z^"^"^^ f . 
Hence, in the case of the torus, jqcp is 



^QCP 



3 
2 



= In 6. 



(3.21) 



2. The compactified boson wave function 

We can also use this approach to compute the entangle- 
ment entropy for the compactified boson wave function 
(the quantum Lifshitz state) discussed in the previous 
Section. However, unlike the explicit computation of the 
boson determinant presented in the previous section, a 
computation that can be done for any compactification 
radius i?, the boundary CFT approach we are using in 



this section only applies for a rational CFT. This restricts 
the compactification radius to be such that is a ra- 
tional number. (The general case can be regarded as a 
limit.) 

It is now straightforward to compute the entangle- 
ment entropy using Ea. (|3.1ip . For this case we find 
IQCP — ~ In 5*0 — Ini?, consistent with the results of 
the preceding section. 



3. Quantum loop models 

Quantum loop models are two-dimensional quantum 
systems whose Hilbert space is spanned by states la- 
belled by loop configurations (or coverings) of a two- 
dimensional lattice. We will denote by {£} the set of 
these configurations. Conventionally, this set of states 
are taken to be a basis of the loop Hilbert space, and 
hence they are assumed to be linearly independent, com- 
plete and orthonormal, (with respect to the naively de- 
fined inner product.) 

Quantum loop models were originally proposed 
as candidates for time-reversal invariant topological 
phases Wave functions in the Hilbert space of 
(multi) loop configurations have the form 



l*(x,<i)) 



(3.22) 



Here N[C] is the number of loops in state (configuration) 
£, L[C] is the length of loop in the configuration, d is the 
"loop fugacity" , and x is the weight (fugacity) of a unit 
length of loop. 

The candidate wave functions of a quantum loop model 
in a putative topological phase depends on the loop con- 
figuration but not on the length of the loops. The sim- 
plest such state is the "d-isotopy" (multi) loop wave 
function"^!^ 



(3.23) 



obtained from \^[x.d)) by setting the fugacity of the unit 
length of loop x — \. This is a generalization of Kitaev's 
"Toric Code" wave function^^ [d — \), i.e. a. 1^2 gauge 
theory deep in its deconfined phase in 2 -I- 1 dimensions. 
Another limit of interest is the "fully packed" state 



oo,rf) / 



hm yx^Md^[^l|£) 



(3.24) 



obtained by setting x — s- cx), which forces the constraint 
that the loops cover the maximal allowable set of links 
on the lattice. 

With the naively defined inner product, the norm 
squared of the d-isotopy state Eq. (|3.23p . is 



<2N[C\ 



(3.25) 
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which is the same as the partition function of a 2D clas- 
sical loop model on the same lattice, with a weight d? 
per loop. Likewise, the norm squared of the fully packed 
loop state |^'(oo,d)) is the partition function Z{oo,(P') of 
the classical fully packed loop model, with fugacity c?^, 
on the same lattice. 

The partition functions of classical loop models on a 
2D lattice have been studied extensively, particularly 
on the honeycomb lattice (for a detailed review see 
Refs . j37ll6^l63| . ) In the fully packed limit, the partition 
function Z{oOt(P) is critical for d < \f2. The universal- 
ity classes of the fully packed loop models (on the hon- 
eycomb lattice) are rational unitary CFTs only for d — \ 
(the S'[/(2)i RCFT) and d = ^2 (the S'C/(3)i RCFT). 
For finite x, the partition function for the dense loop gas 
Z(x,c?^) is also critical for d < \/2. The universality 
classes are again rational unitary CFTs only for d = 1 
and d = \pl. The fixed point for the case d = 1 is equiv- 
alent to the statistics of the proliferated domain walls 
of the classical 2D Ising model at infinite temperaturei^ 
For d — a/2 the dense and dilute loop gases have the same 
critical theory, the Kosterlitz-Thouless critical point, and 
hence also the SU{2)x RCFT. 

We can now use the result in Eqs. (|3.12p and (|3.13p to 
compute the universal term of the entanglement entropy 
for the loop wave functions with d = 1, V2, on a cylinder 
with fixed boundary conditions (for the loops) . The mod- 
ular S'-matrices are known, ^'^''^^'"^^ and the needed 5° ma- 
trix elements are 5'° — for SU{2)i and SU{3)i, 
respectively. The universal term "fQCP of the entangle- 
ment entropy for each case is "fqcp = In \/2, In \/3, — In 2 
for the fully packed state at d = 1 (and also for the loop 
gas at d = \/2), the fully packed loop state at d = \/2, 
and the dense loop gas at d = 1 (corresponding to the 
Kitaev state), respectively. Here we have used a recent 
result on the behavior of of the dense loop model by 
Cardy''* who showed (among many other things) that 
for d = 1 the partition function of the dense loop model 
on the cylinder Z = 2. We will see in the discussion sec- 
tion that this negative value, 7 = — In 2, coincides with 
the direct computation of the topological entanglement 
entropy in the Kitaev wave functionJ^i^i^ 



4-. Quantum net models 

Finally, we will briefly discuss the more interesting, but 
less understood problem of the wave functions for quan- 
tum net models^^'^^'^^'^'' . These states were proposed 
as candidates for a time-reversal invariant non-Abelian 
topological phase. The Hilbert space of quantum net 
models is spanned by the coverings of a lattice by con- 
figurations of nets, i.e. branching loops (with trivalent 
vertices). An interesting example is the chromatic poly- 
nomial state<^ In this state, the nets are regarded as a 
configuration of domain walls of a Q-state Potts model. 
The weight of a given state \C) is the chromatic poly- 



nomial xq[>C] of the configuration. The chromatic poly- 
nomial counts the number of ways of coloring regions of 
the lattice separated by domain walls of a Q-state 2D 
Potts model. They were first introduced in the computa- 
tion of the low temperature expansion for the 2D Potts 
models (see, for instance, Ref. j42|.) For non-integer Q, 
the chromatic polynomial can be computed by an itera- 
tive procedure. The 2D Potts model is known to have 
a critical point for Q < 4. 

Following Ref. [s^ , we consider the norm of the chro- 
matic polynomial state with Q < 4. In order to compute 
the norm, we have to square the weight, resulting in a 
partition function involving the sum of the square of the 
chromatic polynomial. It is then natural to ask for a 
value of Q such that x%[^] °c XQchI'^Ij fo'' some Qcff- 
Then the nets will be critical provided Qcff ^ 4. It turns 
out^^ that, up to a suitably chosen fugacity for trivalent 
vertices^, this property holds only for ^/Q = ^ the 

Golden Ratio, with Qcff = 2 + < 4. Thus, for this 

state the nets are critical. 

This case is interesting for several reasons. One is 
that strong arguments^ suggest that it is possible to 
define for this wave function an excitation (a defect) 
which is denoted by r, a Fibonacci anyon (not to be 
confused with the modular parameter!) with the fusion 
rule, TXT = 1+T. Fibonacci anyons are of prime interest 
in the topological approach to quantum computation;^ 
However, for this approach to work it is necessary that 
this state should describe a topological state, which re- 
quires that its local excitations (not the nets) be gapped. 
Fendley^ has recently given strong arguments that im- 
ply that this state, with the naive inner product we use 
here, is not topological but a quantum critical state. 

Another feature that makes this state interesting is 
that the correlations encoded in the norm of the state 
for y/Q = are described by a RCFT, the mini- 

mal model of the Friedan-Qiu-Shenker— series of unitary 
RCFTs at level m = 9, with central charge c = This 
minimal model has a large number of primaries (36) and 
has not been studied in detail. Nevertheless, its modu- 
lar S"- matrix is known (as it is for the entire series'^A). 
Although to the best of our knowledge the boundary 
CFT of this minimal model has not been investigated, 
we conjecture that the boundary state corresponding to 
the fixed boundary condition is the analog of the state 
|0) in the 2D critical Ising model (the m = 3 member of 
the same series.), i.e. the state in the conformal block 
of the identity.'^^ Thus, if we consider this state on a 
cylinder with fixed boundary conditions, the entangle- 
ment entropy for observing only half of the system, has a 
universal term "/qcp of the form given in Eq. (j3.13[) . and 
hence is given in terms of the Sq element of the modular 
5-matrix of this RCFT:^! 
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IV. CONCLUSIONS AND DISCUSSION 

We have shown that at 2D conformal QCPs (with dy- 
namical exponent z = 2), the entanglement entropy for a 
region with a smooth boundary quite generally has uni- 
versal finite contributions which we denoted by "/qcp- 

Sqcp = + 7QCP- 

We studied the universal nature of Jqc p with two com- 
plementary approaches for large classes of 2D conformal 
QCPs: First for the QLM universality class, we calcu- 
lated 7QCP explicitly in terms of the partition functions 
(that of compactified boson) associated with the norm 
squared of the wave function. Later we used known re- 
sults from boundary CFT to show that Jqcp is deter- 
mined by the detailed structure of the associated RCFT 
encoded in the modular S'-matrix and the OPE fusion 
coefficients for the primary fields. We also applied this 
general results to compute jqcp in several systems of 
interest: the quantum Lifshitz model, the generalized 
quantum dimer and quantum eight-vertex models, and 
quantum loop and net models. 

However, we showed (c.f. Ea. (|3.13|) ) that for a general 
conformal quantum critical point, whose ground state 
wave function is given by the Gibbs weights of a Eu- 
clidean rational unitary CFT, the universal term jqcp 
is determined by the modular S'-matrix associated with 
the norm squared of the wave function. Thus, the mod- 
ular S'-matrix of the topological phase and that of the 
wave functions of 2D conformal quantum critical points 
have a conceptually different origin. 

We note that while our result for the entanglement en- 
tropy has the same form as the entanglement entropy for 
a topological phaser^'-^ the finite universal terms jqcp 
and 7topo have a different origin and structure. In the 
case of a topological phase, 7topo is in general determined 
by the modular S-matrix of the topological field theory 
of the topological phase, and it is given in terms of topo- 
logical invariants of the effective topological field theory 
that describes this phase J^iiSii^ This modular S-matrix 
governs the transformation properties of the ground state 
within the degenerate ground state Hilbert space of the 
topological phase under modular transformations on a 
torus: T — > — 1/t, where r is the modular parameter of 
the torus^i. On the other hand, for 2D conformal QCPs 
whose ground state wave function is given by the Gibbs 
weights of a Euclidean rational unitary CFT, the univer- 
sal term jqcp is determined by the modular S-matrix 
associated with the norm squared of the wave function 
and the S'-matrix connects between different boundary 
conditions. Hence the roles of the modular S-matrix 
in the computation of the universal C(l) terms to the 
entanglement entropy have conceptually different origin. 
Moreover, Jqcp and ^topo enter with opposite signs in 
their contributions to their respective entanglement en- 
tropies. In fact, in all the cases we looked at we found 
that jQCP > 0, except for the Kitaev state which is topo- 



logical, and we recovered the known result. (It is unclear 
to us how general this difference actually is and, more 
importantly, if it has a deeper meaning.) In any case, 
the fact that the entanglement entropy has the universal 
form of Eq. (|1.2[) has led to the widespread assumption 
that this scaling is a signature of a topological phase. 
However we have shown here that this is not necessar- 
ily the case as this scaling is also obeyed at conformal 
quantum critical points in 2D. 

It is also interesting to note the striking similarity of 
the structure of Eq. (|3.12p (with its dependence on the S- 
matrix and the fusion rules) with the results of Fendley, 
Fisher and Nayak^^ for the change in the entanglement 
entropy of a 2D topological fluid, a fractional quantum 
Hall state, by the action of a point contact. Recently, 
Refs. [69ll70 | found finite universal terms in the entangle- 
ment entropy for 1 + 1 dimensional CFTs with a similar 
structure to what we found here in 2D conformal QCPs. 
Calculations of quantum fidelity in ID also find a simi- 
lar structure F^i^ Recent work by Li and Haldan&S also 
raises the interesting possibility of computing the entan- 
glement spectrum for a theory with a wave function de- 
scribed by a known CFT, but this is beyond the scope of 
this paper. 

Finally, given the close connection between the uni- 
versal piece of the entanglement entropy jqcp and the 
Affleck-Ludwig entropy of the associated 2D classical par- 
tition functions it is interesting to inquire if "fgcp may 
flow under some perturbation. Clearly this cannot hap- 
pen under the action of a boundary perturbation (as in the 
Afflcck-Ludwig case) as that would require one to make 
a physical change of the wave function on the boundary 
r, rather than a measurement. However, it is interesting 
to consider instead how the entanglement entropy (and 
in particular the finite term jqcp) would evolve as one 
perturbed the (bulk) system either by a finite non-zero 
temperature into the quantum critical regime, or by a 
relevant operator that drives the system into a nearby 
topologically ordered phase that can be accessed by lo- 
cal perturbationsi^i2£i22i^i^ and to investigate possible 
connections with RCFT,2ii2^i2^ 
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